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Real-time nonequilibrium quantum dynamics of electrons in double-dot Aharonov-Bohm (AB) 
interferometers is studied using an exact solution of the master equation. The building of the 
coherence between the two electronic paths shows up via the time-dependent amplitude of the AB 
oscillations in the transient transport current, and can be enhanced by varying the applied bias on 
the leads, the on-site energy difference between the dots and the asymmetry of the coupling of the 
dots to the leads. The transient oscillations of the transport current do not obey phase rigidity. 
The circulating current has an anti-symmetric AB oscillation in the flux. The non-degeneracy of 
the on-site energies and the finite bias cause the occupation in each dot to have an arbitrary flux 
dependence as the coupling asymmetry is varied. 

PACS numbers: 73.23.-b, 73.63.-b 



I. INTRODUCTION 

Coherence of electronic transport through mesoscopic 
junctions has been studied intensively in nanoelectronic 
systems. In particular, the interference of electron waves 
has been visualized in Aharonov-Bohm (AB) interferom- 
eters via the AB oscillations of the conductance of a ring 
placed between two leads. Following the electron injec- 
tion from the leads into the ring, the electrons undergo a 
nonequilibrium transport process before a steady inter- 
ference pattern is reached. While the steady-state AB 
interference has been largely explored in the literature, 
the real-time dynamics of electronic transport in AB in- 
terferometers has not yet been fully understood. In this 
paper we study this dynamics in double-quantum-dot AB 
interferometers addressing the transient AB interference 
under various tunable parameters of the system. 

The study of the wavy nature of electronic transmis- 
sions has been mainly focused on the complex amplitudes 
of the transmitted electrons in the scattering approach*^ 
The archetype model contains a single quantum dot sit- 
ting on one of the two arms of the AB ring. A quan- 
tum point contact (QPC) placed nearby the quantum 
dot has been used to study the effect of a which-path 
detection^ A single-dot AB interferometer has been re- 
alized in a closed geometry^ and also in an open one£ 
Phase rigidity in a two-terminal geometry has been ex- 
perimentally discovered^ 6 , and theoretically explained/I 
The effect of electron-electron interactions on transport 
through AB interferometers has also been explored^ A 
review on the early progress can be found in Ref. [Toj . 
Extracting the transmission phase from the AB oscilla- 
tions is another main issue. The way continuous phase 
shifts (as opposed to phase rigidity) of the AB oscillations 
can be induced by breaking the unitarity of the scatter- 



ing matrix and the way such phase shifts depend on the 
properties of electron losses have been investigated] 11 ' 12 
Likewise, ways of extracting both the amplitude and the 
phase of the intrinsic transmission amplitude from the 
measured conductance without opening the interferom- 
eter have been suggested^ The studies of AB inter- 
ferometers with two quantum dots placed on the two 
arms of the ring have been focused on different issues, 
such as the flux- dependent level attraction^ the effect 
of intradot Coulomb interactions ) 15 ' 16 inelastic scattering 
with phonons,— as well as extracting transmission phases 
from the current measurements using QPC placed next 
to the one of the quantum dots^ 

The above investigations are concentrated mainly on 
steady-state properties of quantum-dot AB interferom- 
eters. In this paper we consider the transient transport 
behaviors in this system. We consider a double-quantum- 
dot AB interferometer as sketched in Fig. [TJ where a 
single active charge state on each dot is assumed and 
electron-electron interactions are ignored. In a recent 
work^ some of us have studied the electron dynam- 
ics in this system under the condition of identical on- 
site energies of the dots and symmetric couplings to the 
leads. In that study, a phase localization phenomenon 
has been found. In the present paper, we systemati- 
cally explore the general transient transport dynamics 
with non-identical on-site energies on the dots and asym- 
metric couplings to the leads. In particular, besides the 
search for the dynamical flux dependences of the tran- 
sient net current, we also examine the flux dependence 
of the transient electronic occupation in each dot and the 
transient circulating current. The electronic occupation 
on each dot can be measured in experiments and contains 
rich information about the transport processes. The rel- 
atively large circulating current at zero or small bias may 
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provide new insights into electron coherence during the 
transport. 
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FIG. 1: A schematic sketch of the system: the Aharonov- 
Bohm interferometer, consisting of two single-level dots, is 
connected to a source and a drain set at different chemical 
potentials, p L and p R , respectively. The interferometer is 
threaded by a magnetic flux $ measured in units of the flux 
quantum $o = hc/e. 



Here is a summary of the main results we obtain. By 
setting the two electronic leads (as the reservoirs) at ther- 
mal equilibrium initially with no excess electrons on the 
double dot, we monitor the time evolutions of the elec- 
tronic charge occupation, the transport net current and 
the circulating current. When the two on-site energies on 
the dots are identical (namely the double dot is degener- 
ate), regardless of the coupling asymmetry to the leads, 
we find that the total electronic occupation on the double 
dot and the net current are always symmetric in the flux, 
while the occupation difference between the two dots and 
the circulating current are always anti-symmetric in it. 
We also find that the times needed for the total occu- 
pation to reach its steady-state values are much longer 
near zero flux, compared with the case where the flux 
value is away from zero. By breaking the degeneracy 
of the double dot, the net current is allowed to break 
phase rigidity transiently at any bias. The flux depen- 
dence of the total occupation number changes arbitrarily 
as the the coupling asymmetry is varied at finite biases. 
The non-degenerate double dot coupled asymmetrically 
to the leads also drives the circulating current slightly 
away from an anti-symmetric flux dependence immedi- 
ately after the current is switched on, but it then quickly 
becomes completely anti-symmetric in the flux. 

The rest of the paper is organized as follows. In Sec. II 
we outline the basic formalism describing the noncqui- 
librium electronic dynamics for nanoelectronic devices 
in general, and for the double-dot AB interferometer in 
particular. In Sec. Ill we present analytical expressions 
for the electronic occupations, the transient net current 
and the circulating current. In Sec. IV we consider the 
steady-state limit, reproduce known results for this sys- 
tem and compare them with ours. In Sec. V we numer- 
ically demonstrate the transient flux dependence of the 
electronic occupations and currents. Finally, conclusions 
are drawn in Sec. VI. 



II. BASIC FORMALISM 

In this section we give a brief introduction to the 
nonequilibrium quantum theory that can describe tran- 
sient quantum transport and quantum coherence in na- 
noelectronic system o 20 ' 21 and then apply it to the doublc- 
quantum-dot AB interferometer considered in this paper. 

The Hamiltonian of the prototypical nanoelectronic 
system we consider can be written as 



n = n s + He + n T 
,t. 
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where ~H S = y,- Eijaja^ is the Hamiltonian of the central 
system with i, j labeling the electronic levels in the dots, 

H E = J2ak e ak c \k c ak ls tne l ea d Hamiltonian with a la- 
beling the leads and k denoting the states in the leads, 
and Ht = J^iakl^akC^a,. + H - c -] describing the tunnel- 
ing between the dots and the leads. Here a\ (aj and 
c ak (c a k) are the electron creation (annihilation) opera- 
tors for electronic levels i and k in the dots and in lead a, 
respectively. En = E t is the energy of level i, Eij (i / j) 
is the tunneling amplitude between the different levels 
in the dots, and Vi a k is the tunneling amplitude between 
the dots and the leads. Electron-electron interactions are 
ignored. 

Since the central system is open to the electron reser- 
voirs (via the leads) , its nonequilibrium dynamics is nat- 
urally described by the reduced density matrix p(t) which 
is denned by tracing over the states of the leads, 

P (t) = tr E p tot (i) = tr E [ e -' m ( t - to )p tot (< )e m(t - to) ], (2) 

where ptot(i) is the total density matrix of the central 
system plus the leads. Electronic occupations on the dis- 
crete electronic states in the dots can be read from p(t). 
The electronic transport through the central system is 
characterized by the currents flowing from the leads into 
the dots, defined by 7 a =-e^ Y^kea ^tot[c ] ak CakPtot{t)] 
for lead a. This can be further decomposed into sepa- 
rate contributions through each dot: 



ha = ie ^ tr tot[Via k cl k aiptot{t) - H.c.]. 



(3) 

In Eqs. (|2|) and (j3)), trg and tr to t denote the traces over 
the states of the leads and the total system, respectively. 
Throughout the paper, we use units in which h= 1. 

As usual, we assume^ 2 - that the central dot system is 
initially decoupled from the leads, and the leads are ini- 
tially at thermal equilibrium with the chemical poten- 
tial p a and inverse temperature f3 = 1/fcgT for lead a, 
whose Fermi distribution function is given by / Q (e) 

• + 1]. Then the exact equations governing 

the time evolution of the reduced density matrix and the 
transient currents ar o 20 i 21 



dt 



p{t) = -i[n s ,p(t)] (t) +£- a (t)]p(t), (4a) 



ha{t) = e tr a [£ t + Q (i)] = -e tr s [£" (t)p(t)}, 



(4b) 
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where the superoperators Cf a {t) are expressed explicitly 
by 

££(*)p(*) = + 4p(*K] 

+ «mj(i)a|aj/o(i) + H.c.j, 

^ (*)/»(«) = E { a «jW [M*)<4 + 

+ K mj (t)ajp(t)a\ + H.c.j, (5) 

and tr s is the trace over the states of the dots. The first 
term on the right hand side of Eq. (|4a[) is the renormal- 
ized Liouvillc operator of the central dot system. The 
second and the third terms, expressed in terms of the 
superoperators, are non-unitary. The non-unitarity is in- 
duced by electronic dissipation and fluctuation processes 
due to the couplings of the central dot system to the 
electronic reservoirs. The transient transport current is 
determined from the non-unitary dynamics, as shown by 
Eq. (jib)) . Equations (|4|) and ([5]) form the basis of the 
nonequilibrium description of quantum coherence and 
quantum transport in mesoscopic systems. 

The time-dependent dissipation and fluctuation coef- 
ficients in Eqs. (O, K a (t) and \ a (t), are explicitly de- 
termined by the nonequilibrium retarded and correlation 
Green functions of the dot system, denoted here by u(t) 
and v(t)r^ via the relations 



K a (t) = / dTg a (t,T)u(r)u 1 (t), 



(6a) 



X a (t) = / dr {g a (t,T)v(r) - g a {t,T)u(r)} - K a (t)v(t). 

(6b) 

The nonequilibrium retarded and correlation Green func- 
tions of the dot system obey the following dissipation- 
fluctuation integrodiffcrcntial equations of motion 



d_ 

df 
d_ 
di 



u(t)+{Eu(t)+ { dr'g(T - t')u(t') = , (7a) 

J to 

v(t) + iEv (t)+ / dr 1 ' g(r - t')v(t') 
J t 

= f dr'g(T - t')u{t') , (7b) 



subject to the conditions u(to) = I,v(to) = with to < 
t < t , and u(t) = u^(t — r + i ) is the advanced Green 
function. Here E is the on-site energy matrix of the dot 
system, g = J2 a 3a and g = J2 a 9a are the self-energy 
corrections due to the coupling to the leads: 



fl«(r) = J ^r a (u)e-*", 
9a(r) = / ^/ Q ( W )r Q ( W )e— . 



The spectral density T aij (u) = 2irJ2 kea V *ak V j<*kS(u - 
e ak ) summarizes all the non-Markovian memory effects 
of the electron reservoirs on the dot system. 

The correlation Green function v(t), Eq. (|7bp . has a 
general solution in terms of the retarded Green function 
u(t), 

v{t) = / dn / dr 2 u(T - n + t )g{n - t 2 )u(t 2 ). 
Jto Jt 

(9) 

From the master equation, Eq. (|4"aj) . it is easy to find the 
single-particle reduced density matrix in terms of u(t) 
and v(t): 

P$(t) = tr[a]a lP (t)} = u ii ,(t)p^ l (to)u] IJ (t) + ««(*), 

(10) 

where p[]\to) is the initial single-particle reduced density 



(8b) 



matrix of the dots. The currents Eq. (|4b[) can then be 
explicitly expressed as2i 

ha{t) = -2eRc / dr{g a (t - t)v(t) -g a (t- t)u(t) 

Jto 

+fla(t-r)ti(r)pM(to)« t (*)}«■ 

(11) 

This expression is consistent with the result obtained 
from the Keldysh Green function technique, except 
that the initial state dependence [the third term in 
Eq. pop ] is usually ignored in most of the Green func- 
tion treatments^ (see the explicit derivation given in 
Ref. (H|). If the dot system is initially empty, namely, 
Pij\to) — 0, the transient electronic occupations and cur- 
rents can be further simplified: 

p$(t)=v ij (t), (12a) 
Iia(t) = -2eRe / dr{g a (t - t)v(t) -g a (t- t)u(t)} ... 

Jto 

(12b) 

Thus, solving Eq. ([7a| and using Eq. ([9]), we can obtain 
the full information of the transient quantum transport 
dynamics. 

To be specific, we consider in this paper a double- 
quantum-dot AB interferometer schematically plotted 
in Fig. [TJ where each of the quantum dots has a sin- 
gle active electronic state. Then the energy matrix E 
in Eq. ((7|) becomes a 2 x 2 matrix. We also do not 
consider the inter-dot tunnel coupling, namely, E\ 2 = 
E 2 i = 0. The AB magnetic flux is embedded in the 
tunneling amplitudes between the leads and the dots: 
VjLk = V jLk e- l ^ L and V jRk = Vjane**** with the rela- 
tion (j)iL—<t>2L+<t>XR—<i>2R = <t> = 27T<I>/4>o, and <I>o = hc/e 
is the flux quantum. Thus the spectral density involving 
explicitly the threading magnetic flux is given by 

T aij {u>) = 2tt^ V ia kV joik e ±i ^"-^S{L0 - e ak ) , 

(13) 
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where the +(— ) sign is for a = L(R). With the above ba- 
sic formulation, we are able to explore the nonequilibrium 
electronic dynamics in this nanoscalc AB interferometer. 



III. ANALYTICAL SOLUTIONS 

We exploit in our calculations the ubiquitously-used 
wide-band approximation, in which the spectral density 
is assumed to be frequency independent. In general, the 
magnetic phase can be characterized by two variables, 
the magnetic flux threading the ring, 4> = 4>l + <j>R, an d 
the difference, i.e. the gauge degree of freedom, x = 
c f >L ^ R , where 4> a = 4>\ a — 4>2 a - Correspondingly, the 



1 



spectral density is reduced to r Q = T a ( 

where the upper (lower) sign is for a = L{R). The time- 
dependent self-energy correction to the retarded Green 
function of the electron in the double dot is given by 



g(r) = 5(t) 



e L 4> 
T 



(14) 



withTj = [rcos(0/2)±^rsin(0/2)]. Here T = 1^+1^ 
and 5T = Tl — Tr characterize the strength and the 
asymmetry of the coupling to the leads, respectively. 

Even for the most general case of a non-degencrate 
double dot asymmetrically coupled to the leads, the solu- 
tion of Eq. (|7a|) can be found analytically (taking to = 0) : 



u(t) = u q (t)(t - u p (t)p{<I>, x) • <t . 



(15) 



Here a = (a +1 a^,a z ) is the vector of the three 
Pauli matrices, and <jq = I (the identity op- 
erator). We have introduced a fiux-dependent 
and gauge-dependent polarization vector p(</>, x) = 



(p_(0,x),P+(^x),P f (0,x))=(ie^r+ le-^T^iSE) 
containing all the information on the gauge depen- 
dence, the flux dependence and the dependence on 
the asymmetry of the couplings, with SE = E\ — E-2 
characterizing the non-degeneracy of the double dot 
on-site energies. Here p(0, x) = p(</>, x) an< ^ 
r> = yjv 2 cos 2 (^/2) + 5T 2 sin 2 (</>/2) - SE 2 which is 
gauge independent. Without loss of generality, we set 
E = El + E ' 2 = 0, as an energy reference. Then the 
functions u (t) and u p (r) in Eq. (|15[) are given by 



uo, p (t) 



±e 



(16) 



with 7 ± = |(r±r ), which are also gauge independent. 
Substituting Eq. (|15|) into Eq. ([9]), we obtain the corre- 
lation Green function 

v(t) = [ ^u(t, u) f a {u)T a u% w) , (17) 

where u(t, uj) = Uo(t, w)cto — u p (t, u))p(<f), x) • <? and 

= ( iw -7^)* _ 1 p(i"+7l")t 



uo lP (t,w) 



l ± e 1 



1 



(18) 



The gauge degree of freedom parameterized by x appears 
explicitly in the off-diagonal matrix elements of the re- 
tarded and the correlation Green functions u(t) and v(t). 
However, the physical observables, calculated from u{t) 
and v(t), do not depend on x, ensuring the gauge invari- 
ance of our calculations. 

Explicitly, the electronic occupation on each dot is 
given by the diagonal matrix element of v(t), see 
Eq. (|12a[) . The total occupation number N(t) = ni(t) + 
ri2(t), where rii(t) = vu(t), can be expressed explicitly as 



\u (t,iu)\ 2 + (Tl + 25E 2 )\ 



U p (t,Uj) 



2n' 



SE 



f-{oj)6T^ \u Q {t,u)\ 2 + [T^ + 25E'-{r'-Sr') — sin 



-2{Tl + SE 2 )Re 



Uo(t,u})u p {t,uj) 



2rRc 



ul(t,L0)u p (t,Uj) 



(19) 



and the occupation difference between the two dots, 
5n{t) = ni(t) — ri2(t), is given by 

5n(t)=J |lmK( t)W )^M]{r^/ + ( w ) 

+ [ Jr ^-^(r 2 -,5r 2 )]/„H}. (20) 

Becef ± (u) = f L (u)±fn(u). 



On the other hand, the current passing from the left 
lead to the right one through dot i is given by U = In — 
Iiii. Summing up the two currents through the two dots, 
we obtain the transport net current / = \{Jl — Ir)- 
Combining the current I\ flowing from the left to the 
right through the first dot with the current —I2 flowing 
from the right to the left through the second dot gives the 
circulating current I c = I\ — 1%. Explicitly, the transient 
net current is given by 



5 



I(t) 



?^f+(w){ STRe(u (t,io) - ^-up(t,w)) - 2|u (t,w)| 2 + (r 2 [cos 2 | cos , 



r 

r 



r^r 
T 



sin 2 «i £E 



6r 



sm ( 



£r 2 [sin 2 ^cosc 



Jf 1 



r| + ^ 2 )|^| 



T 2 ST(cos 2 | + 1) + <5r 3 sin 2 | + — / F % j 



2 



u* (t,uj)u p (t,uj) 



g/_( w) {Ra(r«o( t>W )-3±^ 



-u p (t,w)) - 



ff 2 (cos 2 1 + 1) +r 2 sin 2 | 



T^2 2 f <m2 • 2 r r 77.2 

I cos ol sin ois 

2 2 



<5r 2 cos 2 



M*,w)| 5 



1^2 2 r r r 2 • 2 r 1 . r 77.2 
1 COS (H Sill h O-fc 

2 2 



r 2 sin 2 



5r 2 (r 2 + r 2 sin 2 + 25s 2 ) ^ w) 



I 2 + 2r<5r 2 Rc 



Uo(t, U))u p (t, b)) 



(21) 



where the dependencies on x m the Green functions u 
and v are exactly canceled by those of the self-energy 
corrections g and g [see Eq. (|12bj) ]. leaving the current 
gauge independent. The transient circulating current is 
given by 



Ic(t) 



U p (t,Uj) 



(r 2 - sr 2 



+ 5T5E 



.u p (t,uj) 



u p {t, oj) 



dh) . , _ [ ,u p (t,u>) ST 2 



2tt 



■ (r 2 -^r 2 )sin< 



2r ^")-fl^ M )V*'")]l. (22) 



r, 



These dynamical quantities depend on the amount of 
the non-degeneracy <5-E, the coupling asymmetry ST, the 
magnetic flux <j) and also the bias voltage applied on the 
leads through the particle distributions in the two elec- 
tronic reservoirs. The time scales for the transient behav- 
iors of these physical observables are determined by the 
factors 1 /7* = 2/(r±r^) in Eq. (JTSJ), in which the flux as 
well as the coupling asymmetry and the non-degeneracy 
play their important roles. 

The symmetric and degenerate double-dot interferom- 
eter has been widely studied in the literature. This cor- 
responds to SE = ST = 0. Thus r = T| cos f |, and the 
above results can be significantly simplified. Explicitly, 
the total occupation number in the double dot is reduced 
to 

N(t) = r f ^/ + h{m;, w )| 2 + \u p (t,u)\ 2 



21 cos 



(0/2)|Re[u5(t,w)iip(i,w)]}, (23) 



and the occupation difference between the two dots be- 
comes 

Sn(t) =rsm(0/2) , COs( /^, x 
I cos(0/2)| 

^/_(w)Im[«S(*,w)«p(t,w)]. (24) 
The transient net current is simplified to be 

m = /g /_(„){ - r Mt,u)\ 2 -\u p (t,u)\ 2 s . n2 1 

+ VRje[u (t,u}) - |cos(0/2)|up(i,u;)]} 



and the circular current is given by 



I c (t) = -rsin(0/2) 



du> 



cos(0/2) 
I cos(0/2)| J 2tt 



(25) 



/+(w)Im[it p (t,w)]. 

(26) 



In general, the transient flux dependence of the physi- 
cal quantities, Eqs. (|T9)) - (j2"2"j) . for non-degenerate double 
dot with asymmetric couplings to the leads, are neither 
symmetric nor anti-symmetric in the flux. The compli- 
cated flux dependencies are mainly determined by the en- 
ergy splitting SE. When the two quantum dots are set at 
degeneracy, SE = 0, regardless of the coupling asymme- 
try and the finite applied bias, both the total occupation 
number and the net current become symmetric in the 
flux, namely, N(fat) = N(-fat) and I (fat) = I (-fat). 
In contrast, the occupation difference and the circulating 
current become anti-symmetric in the flux: 8n(fa t) = 
—Sn(-fat) and I c (<j),t) = —I c (—fat). However, when 
the degeneracy is lifted, both the symmetric and anti- 
symmetric flux dependencies are transiently present in 
all these physical observables. These complicated flux de- 
pendencies can be simplified by setting an applied bias, 
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lit = eV/2 = —fiR. Under such a bias configuration, 
the terms involving f-(uS) in Eqs. ([20)1 and (j2"2")l van- 
ish since both Im[iip(i, ui)u p (t, uj)) and Im[u p (t, ui)\ are 
odd in uj while f-(u)) is even in uj. The difference in 
the occupations of the two dots then becomes symmet- 
ric in the flux, Sn((/),t) = Sn(—(/),t), and is proportional 
to SE. On the other hand, the circulating current gen- 
erally contains two contributions. One is proportional 
to (T 2 — ST 2 ) sin(</>) and is anti-symmetric in the flux. 
The other contribution is proportional to STSE, and is 
symmetric in the flux. However, we find that the sec- 
ond contribution decays to zero within a time scale of 
a few 1/r. After that time, the circulating current be- 
comes anti-symmetric in the flux, proportional to sin(0). 
Applying the aforementioned bias configuration does not 
affect the existence of both the symmetric and the anti- 
symmetric flux dependent components of the total occu- 
pation and the transient net current. Only in the special 
case of zero bias, the total occupation number becomes 
symmetric in the flux. The transient net current always 
has a non- vanishing anti-symmetric flux dependence for 
arbitrary values of the bias when 8E ^= 0. 

More interestingly, during the transient transport pro- 
cesses, /(</>, t) 7^ I(—(j),t) for the non-degenerate case. In 
other words, it transiently breaks the well-known phase 
rigidity at arbitrary biases. This is easily understood be- 
cause during the nonequilibrium transient processes there 
is no time-reversal symmetry. The time-reversal symme- 
try is the prerequisite for phase rigidity of the linear con- 
ductance of a two-terminal device4&£ Only at steady 
state, as we show in the next section, can Eq. (j2"Tj) re- 
produce this phase rigidity, independent of the value of 
SE. The on-site energy splitting thus plays a crucial role 
for the time-reversal symmetry breaking with respect to 
the flux dependency during the transient dynamics. Note 
also that /(</>, t) ^ I{~<f> 7 1) is a purely transient effect. In 
a steady state, the phase rigidity is preserved. 



IV. COHERENCE AND PHASE RIGIDITY AT 
STEADY STATE 

Before studying the real-time dynamics of electronic 
transport in this double-dot AB interferometer, we de- 
duce the steady-state results from the general transient 
solutions given in Sec. Ill, and compare them with the 
previous steady-state solutions obtained for a degenerate 
double dot coupled symmetrically to the leads. 



A. The general steady-state solution 

Taking the steady-state limit of Eq. ((TO)) , we have 
lim t ^ 

oo u o,p(£i k- 1 ) — 2^ u ~^ Ll) ^ u +( t -* J )) where w±(w) — 
± . . Substituting these solutions into Eqs. (|T9|) - (|22l) . 

we obtain the electronic occupation and transport cur- 
rents at steady state. 



The total electronic occupation at steady state is re- 
duced to 



/OO J 
2 
.oo 27T [uj 2 + ( 7 + ) 2 ] [uj 2 + ( 7 )2] 



dui 
2^ 



/-(«){■ 



+ 7^7 ) 

t)V+( 

5T(uj 2 -^y) 
" 2 + ill) 2 ][co 2 + { 1 -A 2 ] 



5E- 



2 4 



sin(0) 



}• (27) 



[u 2 + ( 1 +) 2 ][uj 2 + ( 1 ;) 2 

At zero bias, only the first term survives. At a finite 
bias, the difference of the particle distributions between 
the two electronic reservoirs can give an additional con- 
tribution to the total occupation when the asymmetric 
coupling or the non-degeneracy are present. Besides, it 
shows that the flux dependence of the total occupation 
has an anti-symmetric flux dependence only when a fi- 
nite bias is applied for a non-degenerate double dot. The 
steady-state solution of the occupation difference is 



5n(<f>) = 



dw' 



TSEf + (uj) + (STSE- 



(r 2 ~sr 2 ) 



sin(c/>) ) /_ (uj 



2tt 



(^ 2 + W) 2 )(^ 2 + (7l) 2 ) 



(28) 



Equation (|28j) shows further that one must have either 
an on-site energy splitting or a nonzero flux under a finite 
bias to generate an occupation difference between the two 
dots. 

The general expression of the steady-state net current 



/ ^1/lM-/hNT(«, 



where the transmission coefficient is 
7"(w,0) 



(T 2 - ST 2 )[uj 2 cos 2 % + (M sinf) 2 ] 



(29) 



(30) 



By taking ST = and 5E = Eq (|30| reproduces 
the results investigated in Refs. |14lil5| . Equation ([2"9")) 
shows that for this two-terminal device the steady-state 
net current (as well as its derivative with respect to the 
bias, i.e. the differential conductance) obey phase rigid- 
ity I(<p) = I (—<f>). The AB flux profile of the steady-state 
net current exhibits a phase shift of only or ir with re- 
spect to all possible variations of the system parameters, 
including the non-degeneracy of the double dot and the 
asymmetry of the couplings to the leads. Correspond- 
ingly, the steady-state circulating current is given by 

/ c (0) = (r 2 -st 2 )x 

duj uj [-£ 8in(0)/+(w) + (SE+§ sin(0))/_H] 



2vr 



UJ* 



hi) 2 } 



(7l) 2 ] 



(31) 
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Both the net current and the circulating current, sec 
Eqs. P? ]) and ([3"T ]) . are proportional to T 2 - ST 2 . Their 
AB oscillation amplitudes decrease upon increasing the 
coupling asymmetry. 



B. Small and large bias limits at zero temperature 

The steady-state occupation numbers and currents, 
Eqs. (|27| -([3T |) . are expressed in terms of integrals over 
the frequency. These integrals can be explicitly car- 
ried out at zero temperature with the bias configuration 
Ht = eV/2 = -fj, R . 

At zero temperature, the total electronic occupation is 
found to be 



1+lST-SE 



sV/2 



STSE T 2 -ST 



- ST- SE-±- 



2 ^ 2 ■ i tan" 
sin r 



eV/2 



nTd 



(32a) 



Assuming a small or a large bias, Eq. (|3"2"j) can be further 
simplified, 



N{4>) -> 1+ 



eV 



5E 



2ar(5_E-(r 2 -<5r 2 )sim 



ST 



SE 2SE5r-(r 2 -ST 2 ) sin 
2r (r a -5r 1! )sm^(0/2)+<5_E 2 



if eV « T 

if e y > r 

(32b) 



when is real, and otherwise 



5n(4>) 




tan 



if eV < T 
if eV > T 



(33c) 

Therefore, when a small bias is applied, the on-site en- 
ergy splitting effectively causes a difference in the occu- 
pations. However, when we apply a bias much larger 
than the energy splitting, the energy splitting becomes 
ineffective in rendering the occupation difference between 
the two dots. The bias setting with respect to the energy 
splitting is thus essential for the control of the occupation 
difference between the two dots. 

Having examined the occupations in these limits, we 
now turn to the currents. The steady-state net current 
at zero temperature is found to be 



(r 2 - ST 2 ) 

■kTT^ 

r 2 <5£ 2 sin 2 § 



cos -- 



47 



tan 



-1 r^7 2 i 



SE 2 sin 2 % 



47; 



1 reV72 n 
2 * tan -1 r 



my 
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(34a) 



For small or large biases, it is further reduced to 



i(<f>) -> (r 2 - st 2 ) 

- .f SE 2 sin 2 I 

X < - - ~ d, , <5B 2 sin 2 t 



(r 2 -5r 2 )sin 2 |+<5_E 2 



if eV <C T 

if e y<r 

(34b) 



Note that since Coulomb interactions have been ignored, 
the screening effect is altogether discarded. By setting 
SE = 0, the total occupation at large bias becomes inde- 
pendent of the flux. In contrast, at small bias the total 
occupation is flux dependent when the coupling to the 
leads becomes asymmetric. 

The occupation difference between the two dots at zero 
temperature reads 



In the small bias limit, the amplitude of the AB oscil- 
lation in the net current increases with the on-site en- 
ergy splitting. However, under a large bias, it shows two 
competing oscillations, cos 2 (</>/2) and sin 2 (0/2), which 
results in a sub-oscillatory pattern over the main oscilla- 
tion of cos 2 ((/)/2), proportional to SE. 

The explicit expression for the steady-state circulating 
current at zero temperature is 



5n((f>) 



SE m (eV/2) 2 + ( 7 ;) 2 



2tiT 
SE 



(eV/2) 2 + ( 7 +) 2 



W*\ 



if is real 
otherwise 

(33a) 



As expected, the occupation difference is proportional to 
SE. The small and large bias limits are 



W) = 



r 2 - st 2 



sin i 



1 

2T <t 

W7\ 



In 



2tt 

(eV/2) 2 +( 7 7) 2 



(eV/2) 2 + (~lJ) 2 

1 ( e y) 2 +r 2 -|ry| 2 
lclIi 2r|r^| 2 



and for small or large biases it reads 



if is real 
otherwise 

(35a) 



6n((f>) 




ifeV^<cr 
if eV > T 



(33b) 



W) 



(r 2 -sr 2 

27iTj, 



sin^ln^r ifeV<.T 
if eV > T 



(35b) 
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when r^, is real, and otherwise 



i c {4>) 



(r 2 -(5r 2 ) 



tan 



-l r^-|r^,| 2 



if eV < T 
if eF > T 



(35c) 



From the above results we find that the circulating cur- 
rent becomes significantly large when the bias is suffi- 
ciently small. In the opposite limit, the large bias drives 
the electron to flow in one direction and the circulating 
motion is then strongly suppressed. 

It is worth noting that in the case of the degener- 
ate double dot at zero flux, the operator A_A-, where 
A_ = -j=(e~ lx / 2 a\ — e lx ^ 2 a 2 ) and \ is the gauge degree 

of freesom, gives a constant of motion>2! H] = 0. 

When one turns on a finite flux, this symmetry is broken, 
and the electronic occupation is changed significantly 
from the value at zero flux. Indeed, if one sets SE = and 
<f> — in Eq. (fT5|) and takes the steady-state limit, one 
obtains N(<p = 0) = 1/2 + |^ tan -1 [^f] at zero temper- 
ature. However, taking the zero flux limit in Eq. (|32[) at 
degeneracy, one finds N{<j> — »• 0) = 1— fp + ^ tan -1 . 
This indicates that 7V(0 = 0) ^ N{<j> -> 0), namely the 
total occupation, changes abruptly across the zero flux 
point. On the other hand, Eq. (|20p shows that at de- 
generacy Sn((f> = 0,t) = lini0_>.o Sn((f), t) = 0, namely the 
occupation difference Sn is continuous across the zero 
flux point. By setting first 5E = and <f> = in Eq. (fT5|) 
and then taking the steady-state limit, compared with 
the limit <f> — > after the steady-state limit is taken, we 
find that both / and I c are continuous as the zero flux 
point is crossed. Thus the abrupt change upon crossing 
zero flux occurs only in the total electronic occupation 
due to the existence of an occupation constant of motion 
at (j) = 0. 

The results presented in this section give the general 
AB flux dependence of the electronic occupation and elec- 
tronic transport in the steady-state limit for the non- 
degenerate double dot coupled asymmetrically to the 
leads. 



REAL-TIME DYNAMICS 



Having the analytical solution for the electronic occu- 
pations and the transport currents in the double-dot AB 
interferometer for an initial empty state, we now examine 
the real-time dynamics of the electrons for various values 
of the on-site energy splitting, the coupling asymmetry 
as well as the externally-applied bias. For simplicity, we 
exploit the bias configuration, fiL = eV/2 = —fJ,R. 



A. Degenerate double dot with asymmetric 
couplings to the leads (5E = but 5T ^ 0) 

When 5E = 0, one can sec from Eq. (|I6|) that 
the time needed to reach the steady-state limit be- 
comes considerably longer as <j) approaches zero. This 
is because the dominant decay factor given by 7^ = 
^(r — T^) in Eq. (|TB]) becomes smaller as = 



r 2 cos 2 (0/2) + (5r 2 sin 2 0/2) > becomes larger, 
when </> approaches zero. However, at <f> = 0, T— T^ = q = 
and the time to reach the steady-state limit, given by 

(7^=0 r 1 = h( T + r 0=o) _1 = T^ 1 * becomes much 
shorter. This is because the applied magnetic flux breaks 
the occupation symmetry associated with the degeneracy, 
as we have discussed at the end of Sec. IV. As a result, 
the total occupation is discontinuous across the zero flux 
point at steady state. In the time domain, this effect is 
manifested as the apparent elongation of the time scale 
for reaching the steady-state limit at small but nonzero 
fluxes. The nonequilibrium occupation dynamics of the 
system with and without a threading magnetic flux be- 
comes therefore significantly different. 



-2 -1 

tJ>/-JX 



1 2 




FIG. 2: Flux dependence of the occupations and currents at 
several different times for SE = 0. The quantum dots are 
initially empty. The difference in the occupation numbers is 
not shown since it remains zero at degeneracy. The black 
short-dashed line is for t = 2/T, the green dash-dotted line is 
for t — 10/r, the blue long-dashed line is for t = 40/r and 
t = 00 is the magenta solid line. The bias is eV = ST, the 
asymmetric coupling is ST = — 0.5r and the temperature is 
ksT — r/20. The parameters used here are also used in other 
figures unless otherwise stated. 

The time-dependent AB oscillations of the total oc- 
cupation, the net current and the circulating current 
are shown in Figs. EJa), (b) and (c), respectively, for 
the degenerate double-dot asymmetrically coupled to 
the leads. The flux dependencies N((j>,t) = N(—<f),t), 
I(4>,t) = I(—(f>,t) as well as I c (t) cx sin(</>) discussed in 
Sec. Ill are shown there. In Fig. Eta), we see that the 
curve for t = 40/r deviates from the curve for t = 00 
for fluxes near zero. A discontinuity of the total occu- 
pation across the zero flux point is shown for the curve 
at t = 00. This demonstrates the occupation symmetry 
breaking by the applied flux at degeneracy as discussed 
in Sec. IV. In Sec. IV, we have also pointed out that both 
the net current and the circulating current are continu- 
ous across the zero flux point. Therefore the long times 
needed for the occupation to reach the steady state near 
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zero flux are not expected for these currents, as shown in 
Fig. 13(b) and (c). 



B. Non-degenerate double dot with symmetric 
coupling to the leads (SE / but <5r = 0) 

We proceed to examine the case with arbitrary on-site 
energy difference on the two dots coupled symmetrically 
to the leads, i.e. SE ^ but ST = 0. 



(al) t=2/r 




(a2) t=40/r 



(a3) t=CQ 



(t>1) 



(b2) 



(b3) 



■■'// \ 



, (c1) 



(C2) 




FIG. 3: Flux dependencies of the occupations and currents at 
several different times with various energy splittings SE. The 
magenta solid lines are for SE = 0, the brown dash-dot-dot 
lines are for SE = 0.15r, the black short-dashed lines are for 
SE = 0.5r and the orange long-dashed lines are for SE — 2F. 



The flux dependence profiles of the occupations in each 
quantum dot, the net current and the circulating current 
are plotted in Fig. [3] for several values of the on-site en- 
ergy splitting. With a small energy splitting, for example, 
5E = 0.15r, we can still see that the times needed for the 
occupations to reach the steady state at fluxes near zero 
are much longer than those at other values of the flux 
(compared Figs. E);a2) and (b2) at t = 40/r with (a3) 
and (b3) at t = oo). However, the occupations at steady 
state are continuous across zero flux when the degeneracy 
of the double dot system is removed, even only slightly, 
as shown by the curve for 8E — 0.15r in Figs.[3la3) and 
(b3). By further increasing SE, the electrons are more 
likely to occupy the level with the lower energy as ex- 
pected, see the curves for SE = 2T in Figs. EJa2), (a3) 
and (b2), (b3). Non-symmetric flux dependencies are 
also observed for the occupations at SE ^ 0, as discussed 
in Section III. 



On the other hand, upon increasing SE to a large value 
(here SE ~ 2r), there is a tt phase jump with respect to 
SE = in the net current, see the orange long-dashed 
lines in Figs.^cl) to (c3). The phase jump upon chang- 
ing SE in the steady-state net current can be easily found 
from Eq. (|2"9"|) . regardless of the value of ST. More inter- 
estingly, the AB oscillation pattern for different SE at 
time t = 2/r, sec Fig. GUcl), has a different phase shift 
as compared to those at a later time, see Figs. |3|c2) and 
(c3). This shows the transiently breaking of the phase 
rigidity in the transient net current, as we have discussed 
in connection with Eq. (|2"Tj) . In contrast, for the case of 
the degenerate double dot, phase rigidity remains at all 
times, as shown in Fig. [2J The behavior of the circu- 
lating current in response to the change of SE is rather 
simple. Its AB oscillations remain proportional to sin</> 
at all times for different SE with ST = 0, sec Figs.[3jdl) 
to (d3), as expected from Eq. 



C. Non-degenerate double dot with the 
asymmetric coupling to the leads (SE ^ and SF ^ 0) 

After examining the effects of the asymmetric coupling 
and the energy splitting separately on the electronic oc- 
cupation and transport dynamics, we next study the ef- 
fects of the asymmetric coupling together with a finite 
energy splitting 



(a2) t=4/r (a3) t=Q3 




FIG. 4: The occupation number of the first dot, the net 
current and the circulating current as functions of flux in 
(al) to (a3), (bl) to (b3) and (cl) to (c3), respectively, with 
SE = 0.5r. The magenta solid lines are for SF — 0, the 
black short-dashed lines are for SF — — 0.5F and the blue 
long-dashed lines are for SF = — 0.8F. 

Figure [4] illustrates the results for the occupations and 
currents under various choices of <5r's with SE = Q.5T. 
The occupation of the second quantum dot is not shown 
for this energy splitting since n-i does not differ much 
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from m. At time t — 2/T, the occupations for different 
asymmetries in the coupling have the same AB oscillation 
phases, see Fig. EJal). At long times, the differences in 
the AB oscillation phases due to different asymmetric 
couplings become more distinct. This demonstrates the 
arbitrary flux dependence of the AB oscillations of the 
total occupation with respect to the change of <5I\ 

FiguresQJbl) to (c3) show that the amplitudes of both 
the transient net current and circulating current decrease 
upon increasing the coupling asymmetry, consistent with 
what we have found from Eq. (f2U)) and Eq. (f3"Tj) . In- 
specting the curves in Fig. 0Jcl), one sees that the flux 
dependence of the circulating current has a small devi- 
ation from the sin </> profile during a short time initially, 
but then it reaches the profile proportional to sin <j>. 

(a-l) t=2/r (gg) t=4/r (a3) t =co 



0.6 





-0.2J . . . 1 J . . . 1 J 

(d1) 102} (d3) 



0.4 



-0.4 




-2-1 1 2-2-1 1 2-2-1 1 2 

<|>/7t §/n §/n 



FIG. 5: The occupation numbers (al)-(b3), the net current 
(cl)-(c3) and the circulating current (dl)-(d3) as functions 
of flux. The magenta solid lines are for eV = 0, the orange 
dash-dash-dot-dot-dot lines are for eV — 0.0125r, the green 
dash-dot lines are for eV = 0.5F, the black short-dashed lines 
are for eV = V and the blue long-dashed lines are for eV = 3I\ 

So far, we have examined the effects of changing the 
asymmetries of the system, namely, SE and ST, at a 
fixed bias. We now turn to study the effects of vary- 
ing the bias. The flux dependence profiles of the elec- 
tronic occupations for different choices of the bias are 
plotted in Figs.[5jal) to (b3), for a given energy splitting 
and asymmetric coupling, SE — 0.5T and ST = — 0.5T. 
When the bias increases, the difference in the occupa- 
tions becomes smaller [see the curves for eV = T and 3r 
in Figs. EJa2), (a3) and (b2), (b3)]. The symmetric flux 
dependence of the net current is transiently broken for 
all biases [see Figs.[3Jcl) and (c2)]. At zero bias, the net 
current goes to zero at steady state, as expected, but a 



finite transient net current is observed [see the curve for 
eV = in Figs.JHIcl) and (c2)]. At finite biases, when the 
net current evolves to its steady-state value, the AB os- 
cillations develop a sub-oscillatory pattern, proportional 
to sin 2 (</>/2), over the main oscillation of cos 2 ((/)/2) [see 
the curve for eV = 0.0125r in Fig. M.c3)]. The anti- 
symmetric flux dependence is maintained for the circu- 
lating current, only the AB oscillation amplitudes vary 
in time for different biases [sec Figs.^dl) and (d3)]. 



VI. CONCLUSIONS 

In this work we have explored the transient quantum 
dynamics of a double-quantum-dot AB interferometer us- 
ing the exact solution of the master equation. We ana- 
lyzed the effects of various tunable parameters of the sys- 
tem, namely, the splitting of the on-site energies on the 
double dot, the asymmetric coupling to the left and the 
right leads and the externally-applied bias, on the time- 
dependent electronic occupations and the net current as 
well as on the circulating current, during the noncquilib- 
rium transient processes. In the steady-state limit, we re- 
cover the results that have been extensively investigated 
in the literature. 

With identical on-site energies on the double dot, re- 
gardless of the coupling asymmetry to the leads, we 
find that the total electronic occupation in the double 
dot and the net current are always symmetric in the 
flux, while the occupation difference between the two 
dots and the circulating current are anti-symmetric in 
it. We also find that the time needed for the total oc- 
cupation to reach its steady-state value is much longer 
near zero flux, compared with the flux values away from 
zero. This is because there exists an occupation sym- 
metry at zero flux, where a discontinuity cross zero flux 
in the total occupation is found. By breaking the de- 
generacy of the double dot, the phase rigidity in the net 
current is broken transiently at an arbitrary bias. By 
varying the non-degeneracy of the double dot and the 
coupling asymmetry to the leads, the total occupation 
has an arbitrary flux dependence at finite biases. The 
non-degenerate double dot with an asymmetric coupling 
to the leads makes the circulating current to slightly de- 
viate from the anti-symmetric flux dependence initially, 
but it then quickly approaches the AB oscillations with 
the fully anti-symmetric flux dependence. The net cur- 
rent shows a sub-oscillatory pattern over the main oscil- 
lation of cos 4> at finite bias. It is also shown that a small 
bias causes a large circulating current whereas the net 
current is negligible. Thus the circulating current may 
provide new insights into electron coherence during the 
transport. 

In short, the splitting of the on-site energies on the 
double dot and the bias configuration applied to the leads 
change significantly the flux dependencies of the transient 
electronic occupations as well as the transient transport 
currents. We hope that experimentally monitoring the 
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transient behaviors will deepen our understanding of the 
electronic dynamics in quantum-dot AB interferometers. 
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